Curvature and torsion in growing actin networks 
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Abstract. Intracellular pathogens such as Listeria monocytogenes and Rickettsia 
rickettsii move within a host cell by polymerizing a comet-tail of actin fibers that 
ultimately pushes the cell forward. This dense network of cross-linked actin polymers 
typically exhibits a striking curvature that causes bacteria to move in gently looping 
paths. Theoretically, tail curvature has been linked to details of motility by considering 
force and torque balances from a finite number of polymerizing filaments. Here we track 
beads coated with a prokaryotic activator of actin polymerization in three dimensions 
to directly quantify the curvature and torsion of bead motility paths. We find that 
bead paths are more likely to have low rather than high curvature at any given time. 
Furthermore, path curvature changes very slowly in time, with an autocorrelation decay 
time of 200 seconds. Paths with a small radius of curvature, therefore, remain so for 
an extended period resulting in loops when confined to two dimensions. When allowed 
to explore a 3D space, path loops are less evident. Finally, we quantify the torsion in 
the bead paths and show that beads do not exhibit a significant left- or right-handed 
bias to their motion in 3D. These results suggest that paths of actin-propelled objects 
may be attributed to slow changes in curvature rather than a fixed torque. 
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Introduction 

The polymerization of cross-linked networks of actin filaments is an important biological 
phenomenon that drives such diverse processes as cell crawling, phagocytosis and 
morphogenesis [H El E] . These complex networks are organized by a number of proteins, 
including nucleation-promoting factors that trigger the formation of nascent filaments 
off the side of existing network filaments [H [21 Hj. Bacteria such as L. monocytogenes 
and R. rickettsii propel themselves during infection by expressing a single protein on 
their cell surfaces, ActA and RickA respectively, which stimulates actin growth by 
activating the branching complex Arp2/3 [5]. Actin polymerization at the cell's rear 
gives rise to a characteristic comet tail that propels the cell through the host cytoplasm. 
Several biophysical models have been developed to describe the connection between 
network growth and the generation of force P El El E] . Despite the extensive theoretical 
and experimental investigation of polymerization forces and velocities [TOl [TT] , less is 
understood about the changes in direction of network growth. 

Comet tails are often highly curved and can propel cells and beads in nearly circular 
or looping paths [121 HSl E]- Phenomenological description of the shapes of the comet 
tails and the trajectories of actin-propelled objects has been used to guage the effect 
of biophysical and biochemical perturbations on the process of network formation and 
maintenance [Ml [IHl [16] . In recent work that examined moving Listeria monocytogenes 
cells confined to two dimensions, Shenoy et al. postulate that a roughly constant torque 
drives cells in looping paths[l7j. A detailed quantification of actin-propelled trajectories 
in three dimensions can yield additional information about all the forces and torques 
that, on average, act on the moving object. Indeed, theoretical treatment of network 
curvature has been used to place limits on the number and placement of pushing 
filaments at a load surface [18]. However, the lack of three-dimensional data on this 
subject has hindered progress in understanding the mechanisms controlling curvature. 

Traditional measurements of object trajectories face two main difficulties. First, 
experiments often track the in-plane objects movements or fluorescence of actin in 
the network, inherently 2D techniques that do not allow one to monitor the full 3D 
complexity of the trajectories. Second, due to the limited depth of field of conventional 
light microscopy, a micron-sized gap between a coverslip and microscope slide is used 
to restrict out-of-plane motion, leading to pseudo-2D observations [HI [5l [131 [19] . 
Motivated by these limitations and an interest in the mechanism governing movement 
direction, we sought to measure the 3D trajectories explored by actin-propelled beads 
in large chambers using three-dimensional laser tracking. 

Methods 

Motility Assays 

Experiments were performed as described previously [20]. Carboxylated polystyrene 
beads 792 ± 23 nm in diameter were coated with the nucleation factor RickA 
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Table 1. Motility parameters for different geometries 
L [/im] V [nm s""*"]* kra/s [/Um^-'-]"!' ( [A^m]''^ t [/xm"-^]^ 

2 58 ± 38 0.22 ± 0.04 11 ± 6 

80 76 ± 10 0.16 ± 0.03 16 ± 9 -0.02 ± 0.03 

* ± ^D; t± SEM 



purified from Rickettsia rickettsii (a gift from Matthew Welch, University of Cahfornia, 
Berkeley). The assay mixture, containing beads in Xenopus laevis egg cytoplasmic 
extract supplemented with an ATP-regenerating mix and Rhodamine-labeled actin, 
was added to a microscope slide immediately after preparation. We used chambers 
of two thicknesses to monitor bead motion in pseudo-2D and -3D environments. 
Approximately 2 fim chambers were created using 2.1 ± 0.1 fim polystyrene beads 
as spacers between a glass slide and coverslip. 80 fim thick chambers were constructed 
by separating a microscope slide and coverslip with double-sided tape. Beads moved 
with a slightly slower velocity, v, in the 2D environment than in 3D (Table [1]), most 
likely due to a viscous interaction with the top and bottom glass surfaces that affect the 
Brownian Ratchet mechanism of actin-based motility [20] . 

Three-dimensional Laser Tracking 

Our instrument uses a single 809-nm wavelength diode laser to measure microsphere 
position using a position-sensitive detector that monitors scattered light in the back 
focal plane [211 ES]. We use a nano-positioning stage capable of moving the sample in 
three-dimensions over hundreds of microns at a bandwidth of approximately 100 Hz. 
In addition, we record fluorescence images concurrent with optical tracking to confirm 
that bead motion corresponds to network growth. 

During an experiment, a feedback routine, implemented in software, allows us to 
track actin-propelled beads over very long distances. 3D photodiode voltage signals 
[22| [23l [23] are sampled at 2 kHz and converted into 3D position in realtime using a 5th 
order 3D polynomial [22]. The measured bead coordinates are then saved to disk and 
used to update the position of the microscope stage such that the microsphere remains 
in the center of the laser focus. Very low laser powers were used to minimize external 
loads that might perturb actin network dynamics. These powers yielded typical forces 
of less than 0.01 fN. 

Data Analysis 

Time derivatives of the x-, y- and z- position versus time records were performing 
using a 4th order Savitsky-Golay filter with a time constant of one second[25]. These 
quantities were then used to calculate the curvature and torsion as defined by equations 
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Figure 1. Example trajectories of 3D bead motion. Data is shown for a bead moving 
in the (a) and SO-fim (b) chambers. The color scale denotes the z-position, the 

position along the microscope optical axis. 



Results and Discussion 

We used a variant of 3D optical-force microscopy to track 0.8-/im diameter beads coated 
with the nucleation-promotion factor RickA moving in Xenopus laevis egg cytoplasmic 
extract. Using 2-/im and 80-/im thick chambers, we created pseudo-2D and -3D 
environments to directly test the effect of the reduction of dimensionality. In the 
2D chambers (n = 13), beads proceeded in characteristically looping paths similar to 
previous observations (Figure[T^) [121 [13]. We observed cases of nearly circular motion as 
well as shapes which more closely resembled a figure-eight. These paths were markedly 
different from those we observed using the 3D chambers {n = 15), where regular looping 
behavior was not seen (Figure [lb). Instead, beads travelled in randomly curved paths, 
rarely resembling a simple circle or helix. 

Using 3D optical tracking, we sought to quantify the shape of the measured bead 
trajectories. Up to an arbitrary translation and rotation, a geometric path in space 
can be fully described by two generalized curvatures, termed the curvature and torsion. 
The curvature, k, measures the deviance of a the trajectory from a straight line, where 
1/k represents the radius of curvature at a given point. Using the coordinates of the 
trajectory with time, r{t), the curvature can be written as 

|r X r| . , 



The torsion, r, measures the deviance of the trajectory from a planar curve and 
can be written as 
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If the torsion is zero, the trajectory hes completely in one plane. Conversely, a positive 
torsion corresponds to a right-handed helicity whereas a path with negative torsion has 
a left-handed helical component. 

We measured the curvature for beads moving in either two or three dimensions 
using chambers of two different thicknesses. The RMS curvature measured for these 
two situations, ~ 0.2/im~^ (Tabled]), are within error of each other and are in good 
agreement with the curvature deduced from previous video-microscopy studies of motion 
of both beads and bacteria of approximately 0.1fim~^ [HI [131 Ej- 

To our knowledge, only one theoretical treatment of actin-driven motion has 
considered the molecular mechanism that produces looping trajectories. Rutenberg and 
Grant [18] showed that a small number of randomly distributed actin filaments pushing 
a load will give rise to an intrinsic curvature due to the non-zero instantaneous torque 
produced by summing the forces from each filament. For a uniformly-coated bead this 
curvature is given by 



where a is the bead radius and n is the number of pushing filaments. In this model, the 
distribution of curvatures depends only on the RMS value, and for a bulk, 3D geometry 
can be written as [18] 



This distribution contains a single broad peak at k = krms/V^ (Figure [2^ dotted 
line), and thus beads are predicted to move with a "preferred" curvature which leads to 
looping paths. 

To directly test this model we calculated the curvature probability distribution 
from our measured 3D trajectories (Figure [2^). Rather than a peaked distribution, we 
find that the curvature distribution decays monotonically, with smaller curvatures, i.e. 
straighter path segments, more probable that curved ones. Our data appears to rule 
out the simple model of Rutenberg and Grant. 

Bead paths most often had low curvatures; fifty percent of the time the radius of 
curvature was greater than 25fim. If this is true, how can we reconsile the appearance 
of looping trajectories such as those shown in Figure [1^ in the 2D chambers? An 
alternative model to that proposed in reference [18] is that while smaller curvatures are 
preferred, network curvature can only change via slow dynamics. If the time to change 
path curvature is greater than the time to go around a loop, tioop = 27i/kv, a path will 
necessarily complete a loop before significantly altering its trajectory. These dynamics 
may be associated, for instance, with the time required to significantly restructure the 
micro-architecture of the filaments. 

We quantified the temporal dynamics of trajectory curvature by calculating the 
autocorrelation of the measured curvatures over time (Figure [2]d). The autocorrelation 
decays exponentially with an average decay constant, (, of 16 ± 9/im for 3D movement. 
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Figure 2. (a) The curvature probability distribution. The probability of measured 
curvatures for all beads using the SO-^m chambers is displayed (bars). The curvature 
has been normalized by its RMS value from each run. For comparison,the distribution 
predicted by Rutenberg and Grant [18^ is shown (red dotted line) . (b) Autocorrelation 
of the path curvature. The autocorrelation of the curvature as a function of path length 
for a single bead trajectory (open circles) and a fit of this data to a single exponential 
decay (red solid line) 



corresponding to an decay time of t decay ~ 210 seconds. In the 2D chambers, the decay 
constant was lower than in 3D, measuring ll±5/im, but corresponding to approximately 
the same decay time, 190 seconds because of the slower measured velocity. 

From this data it is possible to describe why beads appear to move in looping paths 
when confined to 2D. While a path is most likely to be fairly straight at any instance in 
time, there is a small probability it will become highly curved. If a highly curved path 
is initiated, the slow curvature dynamics cause the path to remain highly curved for an 
extended period resulting in loops when confined to two dimensions. When allowed to 
explore a 3D space, where a constant curvature does not lead to planar loops, highly- 
curved path sections merely look "curly" rather than looping. 

What might cause the 200 s time scale, which is much slower than the characteristic 
times for actin monomer addition, ~ 100s~^, or the addition of new filaments to the 
network,~ 10s~^ [6]? While our beads are round and uniformly coated with RickA, 
there may be some residual asymmetry to the spatial distribution of the nucleator. 
Assuming that this density affects the shape and dynamics of network formation, a 
rotation of the bead relative to the network could yield a change in path curvature. 
Therefore, it is possible that this time may correspond to the rotational diffusion time 
for a bead attached to a growing actin network. A second explanation could involve 
architectural changes in the actin network itself. Recent evidence indicates that a 
growing actin network contains memory of past forces it has encountered [11] . The actin 
networks propelling the bead may actively remodel their micro-architecture in response 
to load and other factors. We propose that the 200 s relaxation time we observe in 
the curvature autocorrelation function may represent the time-scale over which actin 
remodeling occurs. 
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Figure 3. The torsion probability distribution. The torsion of measured curvatures 
for aU beads using the SO-^m chambers is displayed. 



In addition to the curvature, we calculated the torsion for beads moving in 3D. The 
distribution of measured torsions was roughly symetric about zero with a small negative 
skew (Figure [3]). The average torsion for all 3D runs was —0.02 ±0.03, not significantly 
different than zero. Small sections of bead movement often appeared helical (Figure [T]d), 
however overall our data shows that beads do not move in consistently helical paths. 
In addition, we measured the autocorrelation of the torsion and found no significant 
correlation over time. The average autocorrelation decayed to zero with a time constant 
equal to the smoothing time constant used to calculate the derivatives of the position 
with time, approximately 5 seconds. 

Our data is markedly different from the right-handed helical tails which have been 
observed with Listeria monocytogenes moving in high-speed brain extracts [26]. These 
trajectories, along with data that moving Listeria cells rotate around their long-axis 
[27] suggest that our RickA-coated beads and Listeria might move by slightly different 
mechanisms. First, Listeria use the nucleation-promoting factor Act A, which might 
produce a slightly different network architecture than that produced by RickA. Electron- 
micrographs of Rickettsia cells show comet tails with filaments arranged in a parallel 
geometry [28] , although beads coated with RickA have been shown to produce dendritic 
networks similar to those generated by Act A 0. Second, Listeria cells have been 
reported to produce a helical surface pattern of Act A [29], which may give rise to a 
network that applies a constant torque to moving cells, producing helical trajectories 
and the looping trajectories seen in 2D [T7]. Our uniformly coated bead would not 
experience this torque. 

Conclusion 

We have presented experimental measurements of the 3D paths explored by actin- 
propelled beads. In a pseudo-2D environment, beads appear to move in looping paths 
as has been shown previously for measurements using video-microscopy. In contrast, 
relaxing the restriction to 2D movement, we find that in 3D beads do not move in 
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looping paths, but still exhibit curved, complex structures. Differing from observations 
of moving Listeria cells, our beads do not move in helical paths in 3D environment, and 
instead we measure an average torsion of zero. 

While this new data places constraints on the physical mechanisms that give rise 
to the movement of RickA-coated beads, it also speaks to the geometric persistence 
inherent in an actin network. The correlation time we observe in the curvature of actin- 
propelled bead paths may be a limiting time for the change in direction of a growing 
actin network. Biochemical events that attempt to change the direction of a growing 
actin network, such as chemotaxis or phagocytosis, must do this in the presence of an 
intrinsic curvature and temporal correlation like that described here. 
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